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Measurements of precipitation pH, SO. and NG\ were taken from 

August 1976 to April 1979 in the Muskoka-Haliburton region of Ontario. Sea level 
geostrophic trajectories were used to follow the air mass associated with the event 
back 48 hours to its 'octant of origin'. Two types of data analysis were carrried 
out, one using the entire data set (composite analysis) and one using only 
precipitation samples consisting of a single event (event analysis). Concentrations 
of H + , SO. and NO.," in event samples from the S, SW, W and NW octants were 
high. Results from the entire data set showed that precipitation events most 
commonly occurred with trajectories from the S and SW octants. Both event and 
composite analyses demonstrated that the majority of H , SO^ and NO., deposited 
at Muskoka-Haliburton was associated with trajectories from the S and SW octants. 



INTRODUCTION 
The acidification of freshwaters by precipitation has occurred in 
Scandinavia (Oden 1976), the northeastern U.S.A. (Likens 1976, Schofield 1976) and 
the greater Sudbury area (Conroy et al. 1975, Beamish 1976). Recent evidence of 
acidic precipitation in south-central Ontario has been given by Dillon et al. (1978), 
who reported mean pH values of 3.95-4.38 at 8 sites in Muskoka-Haliburton. The 



dominant ions in precipitation at this site are H + and SO. " and the deposition of 



H and SO^" " is higher in Muskoka-Haliburton than at any other location on the 
Canadian Shield measured to date (Scheider et al. 1979a). Effects of acidic 
precipitation on the poorly buffered freshwaters of the area are discussed by 
Jeffries et al. (1979) and Scheider et al. (1979b). 

In this paper, we attempt to relate the deposition and concentrations of 
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H + , SO. " and NO," in precipitation collected at Muskoka-Haliburton to pollution 

source regions through the use of air parcel trajectories. The rationale behind the 
study is to investigate the relative contributions from the various source regions in 
order to enable an assessment of proposed strategies aimed at controlling and 
reducing the acidity in Ontario's precipitation. 

Trajectory analysis, to determine the sources of acidity in 
precipitation, has been applied in Sweden (Munn and Rodhe 1971), Norway (For land 
1973) and at sites in the northeastern U.S.A. (Dittenhoefer and Dethier 1976, Miller 
et al . (1978). Several authors have studied airborne SO." in Ontario and its 
relationship to long range transport (Lafleur and Whelpdale 1977, Chung 1978). 
Summers and Whelpdale (1976) related precipitation pH and SO ~ to air mass origin 
in a study performed in August 1976 in eastern Canada. To date, no longer term 
study has been carried out for this area. Here, we discuss the deposition and 
concentration of H + , SO." and NO," in precipitation collected at 5 sites in 
Muskoka-Halburton over the period August 1976 - April 1979 with respect to air 
parcel trajectories. 
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DESCRIPTION OF STUDY AREA 
Precipitation collectors were located in the Muskoka-Halibruton area of 
central Ontario (Figure 1). This area is underlain primarily by Precambrian 
bedrock with a thin, glacially derived overburden. A mixed deciduous/coniferous 
forest covers the area and no large urban or industrial centres exist within 200 
kilometres. 

METHODS 
Precipitation was collected from August 1976 to April 1979 (33 months) 
at five sites using automatic, moisture-activated samplers. The samplers opened 
only during a precipitation event and thus excluded any dry deposition. The periods 
of collector operation are summarized in Figure 2. One sampler (Dorset) was 
electrically heated for winter use and snow samples were thus melted and kept in 
the liquid state. Rainfall depth was measured at each site using a 10 cm diameter 
polycarbonate rain gauge. The depth of snowfall was measured in 1978-79 with a Nipher 
snow gauge. In previous years, snow samples were melted and their volume measured to 
obtain water equivalent depths of snowfall. All precipitation gauges and collectors 

were located in cleared areas about lm above ground (or snow) level. Samples were 
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analyzed for pH and NO.. within two days of collection. Analysis of SO. was 

performed in Toronto within one week of collection. Analytical methods are 

described in Ontario Ministry of the Environment (1975). 

To trace the origin of the air associated with a precipitation event, we 

have used sea level geostrophic trajectories (Rossby and Egner 1955, Munn and 

Rodhe 1971, Forland 1973), computed from maps of the surface pressure field. A 

trajectory segment is derived for each six hourly surface map with the aid of a sea 

l^vel geostrophic wind scale which provides the length and direction of the 

trajectory as a function of isobar spacing and orientation. The air parcel is followed 

back 48 hours in 6 hour increments to its 'octant of origin'. In some cases, the air 

mass associated with an event remained almost stationary for the 48 hours 

preceding the event. Such cases were termed 'local'. 
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Figure 1: Location of. Study Area and Collectors. (Numbers refer to Figure 2) 
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Figure 2: Periods of collector operation. Numbers refer to location on Figure 1. 
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1. Carnarvon 



2. Eagle 



3. Dorset 



4. Red Chalk Lake- 



5. Harp Lake 
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The precipitation collectors are part of a network designed with the 
intent of calculating the bulk and wet deposition of nutrients, ions and heavy metals 
to lakes in the Muskoka-Haliburton area (Scheider et al. 1979 a, Jeffries and Snyder 

in press). For this reason, precipitation samples were collected when sufficient 
volume of sample existed for complete analysis. Depending upon the depth of 
precipitation in an event, one or several events may have contributed to a single 
sample. 

In this study, we carried out two sets of computations - one for the 
entire data set (composite sample analysis) and one for samples consisting of a 
single precipitation event, or several events each associated with the same octant 
(event sample analysis). A summary of the number of measurements within the 
composite and event data sets is given in Table 1. 

(a) Composite Sample Analysis 

The composite data set consists of all data available, including that from 

composite collections in which several events contributed to the sample. We used 

these data to make estimates of the maximum possible deposition (of H + , SO. " and 

NO-." ) which could correspond to each octant. Considering each octant in turn, we 

summed all of the sample loadings associated even partially with that octant. In 

generating maximum possible deposition we were assuming that events from the 
octant of interest alone contributed the total load to the sample. 

An example may serve to clarify the procedure. Table 2A provides a 
data sample for the Red Chalk Lake station. To calculate the maximum possible H + 
loadings for each of the octants we get 1983 + 1085 = 3068 for South, 3315 for South- 
East, 3315 for South-West, 3315 for North-West and 1085 for North-East 

(Table 2B). In the analysis, this procedure was applied to the 
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H , SO- and NO., data for each station. 
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TABLE 1: Number of Measurements in Composite and Event Data Sets 



Station Carnarvon Dorset Eagle Harp lake Red Chalk Lake 



Parameter 



,+ 



H composite 54 56 53 37 40 



SO." composite 



NO.," composite 



47 34 50 38 34 



40 37 40 38 30 



depth-composite 60 58 61 44 43 

H + event 27 24 28 20 18 



SO^' 2 event 29 15 30 21 14 



NO 3 -1 event 21 16 24 20 11 



depth-event 33 26 35 24 18 
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Table 2A: A data sample for the Red Chalk Lake Station (1978) 



Collection Depth Sample No. of Octant(s) H + load 

Period (mm) pH Events of origin (ueq m ) 



Sept. 11-15 24.4 4.09 2 S, S 1983 



Sept. 15-22 43.7 4.12 3 SE,SW,NW 3315 



Sept. 22-29 11.1 4.01 2 NE, S 1085 



Table 2B: Calculation of Maximum Possible H + Loading for various octants 

of trajectory origin (based on data from Table 2A) 

S 1983 + 1085 = 3068 ueq/m 2 

SE 3315 

SW 3315 

NW 3315 

NE 1085 
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As Figure 2 indicates, sampling at the Dorset station generally 

covered the cold months ('winter') while precipitation collectors at the remaining 

sites were made during the warm ('summer') months. We have therefore grouped 

the data into 'summer' (using the arithmetic mean of loadings at Carnarvon, Eagle 
Red Chalk Lake and Harp Lake) and 'winter' (using the Dorset loadings). 

b) Event Sample Analysis 

A second set of computations was based on the 'event' data i.e., those 
derived from precipitation collections consisting of a single precipitation event, or 
several events each associated with the same octant of origin. Initially, we 
computed total event loadings analogous to the composite loadings described 
above. We then computed weighted average concentrations for each octant and for 
each ionic constituent. The weighted average concentration is calculated from: 

E (EL. ./m.) 



w 
a 



h Mi 

E (Ed. ./m. ) 



where: W = weighted average concentration for Octant "a" 
a 

L- = loading from collector j for event i 

m.= number of collectors or gauges operating for the i-th event 

n = number of events corresponding to Octant "a" 

d.. = depth from gauge j for event i 

RESULTS AND DISCUSSION 
The limitations of this analysis are initially stated. They include: 
1. The majority of the precipitation collected is in the form of composite 
samples i.e., samples comprising two or more precipitation events. The 
value of such data to trajectory studies is limited (see, for example, 
Miller et al. 1978). 
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2. The size of the data set for event samples (i.e., those associated with a 
single precipitation event) is small. There are, for example, only one 
'winter' and one 'summer' event sample corresponding to a trajectory from 
the northwest sector. The average concentrations associated with the 
less frequent sectors cannot be considered reliable. 

3. The characterization of acidic precipitation events using a single 
trajectory may be misleading as (a) the storm system involves inflow of 
air from many directions through convergence in the vicinity of the 
cyclonic cell, and (b) the mechanism of acidic precipitation formation 
is a complex amalgam of rainout and washout processes occurring at 
many atmospheric levels possessing different directions of air flow. 
Despite these difficulties, some success has been reported from the use of 
trajectories for tracing source regions (e.g., Forland 1973, Nordo 1975 and 
Wisniewski 1976). 

Table 3 shows the percentage precipitation event samples associated with 
each octant of trajectory origin as compared with the percentage of composite 
samples associated with each octant of origin, for the H + data set. (The SO. ~ and 
NO.," data sets show qualitatively similar distributions). The Table shows that the 
South and South-West trajectory directions are more frequent in the event data set 
while several other trajectory directions are less frequent (particularly East, 

South-East and North-West). This difference in the frequency distribution results 
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in a corresponding difference in the computed loadings of H + , SO. " and NO., - by 

octant of origin. 

The results of the composite sample analysis (Table 4) show the maximum 

possible deposition of H + , SO. ~ and NO.," which could correspond to each octant 

of trajectory origin. In the case of each of the constituents, and both warm and 

cold seasons, the South and South-West octants are associated with the greatest 
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TABLE 3: Percentage of precipitation samples, on which H analysis was done, 

associated with each octant of trajectory origin for event and composite 
data sets 





Summer* 
Event Data (%) Composite Data (%) 


Win 


iter 


Octant of 
Origin 


Event Data (%) 


Composite Data (%) 


NE 


5 






6 


4 


7 


E 









4 


4 


7 


SE 









6 





8 


S 


25 






19 


42 


22 


SW 


46 






29 


34 


24 


w 


14 






13 


k 


13 


NW 


2 






7 


4 


9 


N 


1 






3 


8 


6 


Local 


7 






7 





4 


Number of 
Samples 


18-28 






37-54 


24 


56 



* Mean Frequencies for 4 stations are tabulated. 



-11- 



-2 



-1 



Table 4: Maximum Possible Loadings of H + , SO , and NO, ( 



_2 

imeq m~ ) 



Associated with Octants of Trajectory Origin (based on composite 
sampling data, Aug. 1976-April 1979) 



H + 



SO 



-2 



NO 



-1 





Summer 


• * 


Winter 




(meq/m ) 


(%) 


2 
(meq/m ) 


(%) 


NE 


4 


6 


34 


32 


E 


7 


10 


28 


26 


SE 


7 


10 


16 


15 


S 


26 


37 


42 


40 


sw 


35 


50 


60 


57 


w 


15 


21 


27 


25 


NW 


10 


14 


32 


30 


N 


3 


4 


4 


4 


Local 


7 


10 


13 


12 


Total** 


70 


162 


106 


241 


NE 


6 


9 


14 


29 


E 


5 


7 


14 


29 


SE 


5 


7 


9 


18 


S 


23 


33 


25 


51 


SW 


34 


49 


31 


63 


W 


15 


21 


15 


31 


NW 


8 


11 


12 


24 


N 


3 


4 


0.4 


1 


Local 


10 


14 


4 


8 


Total** 


70 


155 


49 


254 


NE 


2 


7 


9 


30 


E 


2 


7 


9 


30 


SE 


2 


7 


4 


13 


S 


11 


38 


13 


43 


SW 


12 


41 


IS 


60 


W 


6 


21 


7 


23 


NW 


3 


10 


10 


33 


N 


0.4 


1 


1 


3 


Local 


6 


21 


1 


3 


Total** 


29 


153 


30 


238 



*mean of 4 stations 

** The total load is not the sum of the individual octant loadings (see text). 
Total % figures shown are calculated from total loadings shown. 
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maximum possible deposition and the North octant is associated with the lowest 
maximum possible deposition. 

No more than 10% (and likely much less) of the 'summer' deposition of H + , 
SO' 2 and NO " could be associated with any of the N, NE, E and SE octants. 
From the 'winter' analysis only the North octant was associated with less than 10% 
of the deposition of each of the constituents. However, the 'winter' figures are less 
reliable than the 'summer' values because of the greater number of events per 
collection which is reflected in the larger sum of % values (Table 4). 

The results of the event sample analysis (Table 5), although based on a 
smaller data set, give an estimate of the % of total deposition contributed per 
octant of origin. In the 'summer' period, events originating in the S and SW octants 
combined contributed 75% of the H + , 70% of the SO^" and 57% of the N0 3 " 
deposition. 'Local' events contributed 10%, 16% and 27% of the H + , SO^~ and 
NO., respectively. Events orginating in the W octant contributed 9%, 8% and 
10% of the H + , SO' 2 and NO " deposition. Events from all other octants 
contributed 2% or less per octant of the total deposition. In the 'winter' period, 
events from S and SW octants again were most important, depositing 85%, 87% and 
86% of the H + , SO." 2 and NO.,' 1 respectively. The fraction of the total load 
deposited by events from the W octant dropped to less than 5% and events from the 
'local' sector contributed 0%. 

Basically, both data sets show a similar result i.e. that most of the 
deposition of H + , SO ~ and NO ' was probably associated with events originating 
in the S and SW octants. The W and 'local' sectors were also important contributors 
to the total deposition in the 'summer' months. 

Our results agree with those of Summers and Whelpdale (1976) and 
Whelpdale (1978) who reported elevated levels of SO." and H + in precipitation 
falling in southern Ontario to be associated with southerly flows of air. Indirect 
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TABLE 5: Sampled Loadings of H + , SO " and NO ~ (meq m ) 

Associated With Octants of Trajectory Origin (based on event sampling 
data, Aug. 1976 - April 1979) 



Winter 







Summer* 






(meq/m ) 


% 


H + 


NE 


0.8 


2 




E 


0.0 







SE 


0.0 







S 


11.0 


33 




SW 


14.1 


42 




W 


2.9 


9 




NW 


0.6 


2 




N 


0.8 


2 




Local 


3.4 


10 




Total 


33.6 


100 


*v 2 


NE 


0.9 


2 


E 


0.0 







SE 


0.2 


1 




S 


12.5 


31 




SW 


15.6 


39 




W 


3.4 


8 




NW 


0.7 


2 




N 


0.6 


1 




Local 


6.5 


16 




Total 


40.4 


100 


N0 3 _1 


NE 


0.2 


1 


E 


0.0 







SE 


0.1 


1 




S 


4.9 


30 




SW 


4.5 


27 




W 


1.6 


10 




NW 


0.3 


2 




N 


0.3 


2 




Local 


4.6 


27 




Total 


16.5 


100 



♦mean of 4 stations 



(meq/m ) 


% 


1.4 


5 


0.1 





0.0 





12.6 


45 


11.3 


40 


0.8 


3 


1.5 


5 


0.5 


2 


0.0 
28.2 



100 


0.7 


5 


0.0 





0.0 





2.9 


21 


9.0 


66 


0.5 


4 


0.2 


1 


0.4 


3 


0.0 
13.7 



100 


1.1 


13 


0.0 





0.0 





3.5 


40 


4.0 


46 


0.0 





0.1 


1 


0.0 





0.0 
8.7 



100 
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support comes from Lafleur and Whelpdale (1977) and Chung (1978), who concluded 
from Ontario data that elevated levels of SO. in air measured by high-volume 
samplers were associated with S-SW air flows. At Ithaca, N.Y. some 350 km south- 
east of Muskoka-Haliburton, Cogbill and Likens (1974), Dittenhoeffer and Dethier 
(1976) and Miller et al. (1978) all reported that precipitation originating in the SW 
sector was the most acidic. Wolff et al. (1979) found for New York, N.Y. that 
events originating in the W and SW sectors had the lowest pH values. Samson 
(1978) reported that the highest SO. levels in air collected with high-volume 
samplers at three locations in New York were associated with air parcels from the 
SW octant. Many studies point to the highly populated, industrialized S and SW 
sectors as being the major source areas of polluted air, a reflection of the greater 
strength of sources in these directions (see, for example, Clark 1978). 

Table 6 is a tabulation of weighted average concentrations for H , SO^~ , 
and NO./ , based on the event data set. Weighted average concentrations serve to 
characterize the precipitation associated with various source regions. Comments 
on the Table cannot be too extensive due to the relatively small sample size. For 
example, the very high average concentrations for H + in 'winter' corresponding to E 
and NW trajectories are based on isolated events. Much more event data are 
needed before meaningful conclusions can be drawn. It should be noted, however, 
that a significant concentration of each of the 3 ionic species characterizes 
precipitation associated with the S, SW, W and NW directions. Note also that the 
weighted average concentrations in 'summer' are relatively high for 'local' events. 
This circumstance is possibly due to inflow of air from the Sudbury region (200 km 
north-west of Muskoka-Haliburton) associated with slow moving storms. 
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TABLE 6: Weighted Average Concentrations of H + , SO ~ 2 and NO ~ ! 

(ueq/1) (based on event sampling data, Aug. 1976 - April 1979) 







Sumi 


mer* 


Winter 








ueq/1 


number of 
events 


ueq/1 m 


imber of 
events 


H + 


NE 


36 


3 


43 


1 




E 


- 





204 


1 




SE 


- 





- 







S 


81 


11 


89 


10 




SW 


70 


18 


69 


8 




W 


101 


7 


102 


1 




NW 


106 


1 


229 


I 




N 


81 


1 


17 


2 




Local 


83 


4 


- 





*v 2 


NE 


M 


3 


21 


1 


E 


- 





- 







SE 


21 


1 


- 







S 


SI 


10 


38 


4 




SW 


93 


19 


55 


7 




w 


73 


5 


63 


1 




NW 


m 


1 


23 


1 




N 


58 


1 


21 


1 




Local 


138 


k 


- 





N0 3 _1 


NE 


13 


2 


32 


1 


E 


- 





36 


1 




SE 


6 


1 


- 







S 


35 


8 


30 


7 




SW 


31 


15 


48 


6 




w 


45 


5 


- 







NW 


53 


1 


14 


1 




N 


29 


1 


- 







Local 


9 1 


3 


- 


o 



•mean of 4 stations 



-16- 






CONCLUSIONS 
In summary: 
(i) precipitation events most commonly occur with trajectories from the 

South and South-West octants, 
(ii) an analysis of both composite and event data sets for the 33-month period 

shows that the majority of H + , SO^~ and NO-," deposited at Muskoka- 

Haliburton is associated with trajectories from the South and South-West, 

and 
(iii) the average concentrations of H + , SO^~ and NO-," in precipitation 

for the S, SW, W and NW octants are all substantial; precipitation from 

the remaining octants was very infrequent. 
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